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The role of inhomogeneous phases in the phase diagram of the Nambu-Jona-Lasinio (NJL) and 
the quark meson (QM) model is examined. By means of a generalized Ginzburg-Landau (GL) 
expansion it is concluded that the critical point in the mean-field phase diagram of the NJL model 
is in fact a Lifshitz point where homogeneous spontaneously broken, inhomogeneous and restored 
phases meet. This picture is confirmed by a mean-field calculation for inhomogeneous phases 
with a one-dimensional modulation. For the latter it is shown that recent results within lower 
dimensional models can be extended to 3+1 dimensions. Also the respective phase diagram for 
the QM model is presented. 
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1. Introduction and summary 

Until today the phase diagram of quantum chromodynamics (QCD) is subject to intense theoretical 
and experimental investigations (for dedicated reviews see Ref. [|I|]). Since ab initio calculations 
are limited to small net-baryon densities, possible scenarios at moderate temperatures and densities 
are often discussed within NJL-type models \ typically within mean-field approximation. 
Here we limit ourself to results for inhomogeneous ground states in the NJL model ^ and QM 
model [^]. These are characterized by a spatially varying order parameter and have been discussed 
for QCD at least in the large N limit, where they are expected to form the ground state at suf- 
ficiently high densities ^, Related to this they show up in holographic models and in 
the quarkyonic matter picture that suggests a similar structure for QCD. The investigations of 
these phases is however limited, mainly because they are technically much more involved. Within 
the NJL model such phases have been analyzed at vanishing temperatures applying further trun- 
cations as well as for the so-called chiral density wave [|l^, 11 1. In the latter case the order 



parameter is assumed to be a plane wave and it can be solved on mean-field level for vanishing 
current quark masses. Recently, also the effect of small current quark masses has been discussed 



for this ansatz [12]. 



Following Ref. [ ]130 we first address inhomogeneous phases by means of a generalized GL ex- 
pansion, i.e. by an expansion of the thermodynamic potential as an effective action in the order 
parameters as well as in gradients acting on the order parameter. This is a systematic expansion on 
top of a mean-field approximation in the vicinity of a second order phase transition and therefore 
especially suited to explore the region around the critical point (CP). For the NJL model it then 
turns out that the CP is in fact a Lifshitz point, where the first order phase transition in the phase 
diagram of the NJL model is replaced by two second order phase transition lines that border an 
inhomogeneous phase and these two transition lines meet at the CP. 

In order to give a better picture of this finding and in particular to estimate the importance of 
inhomogeneous phases in the phase diagram also away from the CP, we then focus on a complete 
mean-field calculation, following Ref. [[l^], and consider inhomogeneous ground states that form 
lattices of domain-wall solitons. The key observation here is that for the case of one-dimensional 
modulations in the order parameter, the problem can be reduced to an analogues problem in the l-i-l 
dimensional (chiral) Gross-Neveu (GN) model. For this model inhomogeneous phases have been 



investigated for the large N limit Ill5|, [16|, [17|, |18|, |19|], which technically corresponds to a mean-field 



approximation, and basically all ground-states have been classified at least for the chiral limit [pO|]. 
For the case of the GN model it is furthermore possible to introduce finite quark masses and to 



study their effect on the structure of the phase diagram [|16p. The possibility to use solutions from 
lower dimensional models is mainly due to the structure of the mean-field Hamiltonian, which 
is of Dirac-type, and a similar procedure is e.g. not possible in (color-)superconductors. As a 



consequence the investigation of the latter is much more tedious [21]. 

Being able to investigate the role of inhomogeneous phases in the phase diagram of the NJL model 
at least for one-dimensional modulations, we can confirm the picture obtained by the GL expansion 



' We refer to NJL-type models as models that at least in the applied approximation reduce to the NJL model, possibly 
extended by additional point-like interactions, on a technical level. This includes e.g. simplified ansatze for the gluon 
interaction, the use of different regularizations, the instanton liquid model and the quark-meson model. 
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for the vicinity of the CP and the absence of a first order phase transition line in the phase diagram^. 
Furthermore we analyze the relation between inhomogeneous phases and e.g. the strength of the 
first order phase transition (present in the case of homogeneous phases), explore the role of finite 
current quark masses and discuss results in the QM model. The purpose for the latter is twofold: 
On the one hand we would like to extend the analysis to a larger class of models in general, on 
the other hand the regularization of the NJL model for inhomogeneous phases is non-trivial with 
regard to a combined vacuum and QCD phase diagram phenomenology. 



2. Generalized Ginzburg-Landau expansion 



Following Ref. 1 13], we first concentrate on the two-flavor NJL model given by the Lagrangian 

where \j/ is the 4A'y^A'^c-dimensional quark spinor for Nf = 2 flavors and Nc = 3 colors, 7^ are 
Dirac matrices, is the scalar coupling and m the mass matrix for degenerate quarks with current 
quark mass m. For Nf = 2 the matrices x" are the conventional Pauli matrices. In mean-field 
approximation with (v/^^) = — 2^M(x) and (xlriy^z^xif) = 0' the Lagrangian gets replaced by the 
bilinear functional 

^MF = Wr^d, -Mix))Y- (^W^ . (2.2) 

In the case of a periodic condensate with Wigner-Seitz cell V and using the imaginary-time formal- 
ism, we therefore obtain for the mean-field thermodynamic potential as an effective action in the 
order parameter 

n(r,Ai;M(x)) = - Jin f ^Y&Y&W ( I , {^MF+piwfw) 



= -^Trz),c,/y Log [S ^) + y "' ^ ' ^^'^^ 

with inverse propagator 

S-\x,y) = {ifd^-M{x))8^'\x-y). (2.4) 

Since the evaluation of the thermodynamic potential for an arbitrary function is non-trivial due to 
the functional logarithm, we first aim at an expansion in the order parameter M(x). Setting m = 
for simplicity and substracting the leading order corresponding to the thermodynamic potential of 
the unbroken phase, we formally arrive at 

Aa(r,At;M(x)) = -- £ -Trz,,.,/,v (5oM)" + - / . (2.5) 

Here we have introduced the bare propagator Sq = 5|m(x)=o ^ short hand notation for 

TrD,cj,v (SoM)" = f [ ...[ TTD,cj{M{x)So{x,X2)M{x2)...M{Xn)Soix„,x)) . 
Jx Jx2 Jx„ 

(2.6) 

^This picture could of course be modified by the inclusion of color-superconducting phases. 



(M(x) —m 



a 
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The domain of integration for x is [0, 7] x V and [0, j] x for ;t:2, . . . ,x„. In the chiral limit the 
expressions for odd values of n vanish. Furthermore, in order to arrive at a local functional, we can 
expand the condensate around x as 

= £ ^£>"M(x)(x„-x)« (2.7) 

|a|>0 

and can extract the GL functional to any desired order in gradients and order parameter. Neglect- 
ing possible issues with the regularization for the moment, we can go to momentum space using 
SQ{x,y) = TY,„-^^{p^Y^)~^exp{ip{x — y)) , where po = {2n+ l)nT. Treating the magnitude of 
the order parameter and the gradients to be of the same order, it is then a tedious but straightforward 
exercise to work out 

aoL{T,^;M{x)) = ^M(x)2 + ^(M(x)V(VM(x))2) 



M{xy + 5{y M{x)YM{xy + -(am(x)) 

2 



(2.8) 



where 



cx,, = {-lyUMfK-Tj: f + (2.9) 

This expression has to be taken with some caution as the model is non-renormalizable. As a result 
usually a regularization, as part of the phenomenological model, is introduced. Due to this ad hoc 
procedure a generaUzation of the regularization to inhomogeneous phases is often not unique. We 
could therefore take the pragmatic viewpoint that a generalization of any such ad hoc regularization 
procedure to inhomogeneous phases is assumed to be such that total derivative terms (which arise in 
the calculation of the GL functional) vanish. An alternative approach is a regularization scheme that 
does not rely on an homogeneous ground state, e.g. a propertime regularization for the functional 



logarithm in Eq.(|2.3[). In this case it is possible to show that no additional complication due to 
ultraviolet divergencies arise. 

The most interesting feature of the GL functional is that the coefficient of the M(x)'*-term is equal 
to the (VM(x))^-term. As discussed below this will allow for inhomogeneous phases in the regime 
where a4 < 0. It is also worth noting that the GL functional takes a similar form as in the one- 
dimensional GN model [15, [17]. This is not the case in superconductors [22] where the underlying 
dynamics is different, namely coming from particle-particle and hole-hole scattering near the Fermi 
surface instead of particle-hole scattering in presented case. 

With the generalized GL functional at hand we can explore the vicinity near the CP, which is 
defined by = 0:4 = and > 0. Without giving a specific choice of model parameters, we 
assume that such a point in the phase diagram exists and focus on the effect of the gradient terms. 
Specific examples will be subject of the following section. For a4 > we have a conventional 
second order phase transition at a2 = between a dynamically broken {a2 < 0) and a restored 
phase (a2 > 0). Limiting to homogeneous phases this phase transition is rendered first order when 
going through the CP into the regime where a4 < 0. The transition line for this case is given by 
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054 = — 1^160^06/3. In this case, however, also inhomogeneous phases can be expected since the 
(VM(x))^-term is negative and a curvature in the order parameter can therefore lead to a gain in 
free-energy. 



< 



- -(36/5 c(^a^) 

a^- -(8/3 a^a^) 

a4--(4 a^af" 
a^--(16/3 o(,a^)" 
a-, - ; a , > 




0.5 

1/2 2 

sign(a,)(|a2|/ag) [A ] 



Figure 1: Pictorial presentation of the phase diagram in terms of the GL coefficients: The gray domain cor- 
responds to the homogeneous dynamically broken ground state, the shaded gray to the solitonic ground state 
(at least when restricting to one-dimensional modulations in the order parameter), whereas in the transparent 
domain the unbroken phase is preferred. A is an arbitrary scale. Also stated are various lines discussed in 
the text. 



To explore this possibility further we limit ourself to a one-dimensional modulation, i.e. Mix) = 
M(z). The solutions to ^AQ = are actually known from the investigation of one-dimensional 
models [18, 22]. They are expressed (up to an arbitrary shift) in terms of the elliptic Jacobi sn- 
function as 

Mw{z) = Vvqsn{qz,v), (2.10) 

where v € [0,1] and q being a scale related to the maximum of Mid{z) and the extension of a 
soliton in the chosen z direction (both scales are related in our case). For v = 1 we have M\d{z) = 
^tanh(^.x), i.e. a single soliton and for v — > the shape becomes more and more sinusoidal albeit 
the amplitude also goes to zero. From previous investigations it is known that when increasing a2 
from zero we reach a second order phase transition into an inhomogeneous phase with q = Mq and 
V = 1. At this point the free-energy of a single soliton becomes negative leading to its formation. 
By using Mq known from the homogeneous case and checking where ^An|^(x)=MiD(z) changes 

sign at V = 0, we obtain = — \J ^a2(X(, for this point. We arrive at the onset of infinitely far 
separated solitons. Further increasing a2 decreases V until it reaches zero. Since q stays finite 
the overall magnitude of Mid{z) given by ^/vq then vanishes and we find a second order phase 
transition to the unbroken phase. 

In case of a second order phase transition from the inhomogeneous phase to the unbroken phase, 
the value of 04 in terms of 02 0^6 is actually universal also for higher dimensional modulations of 
the order parameter. Since in this case M(x) is parametrically small, we can neglect non-quadratic 
terms in the GL functional. Consequently the variation ^AH leads to a linear partial differen- 
tial equation. We can then optimize the value of 04 by varying the momentum q of the Fourier 



components of M(x) and find 04 = — w |a20J6 for the transition line where |q| = J — j^ 
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The modification of the phase diagram in the vicinity of the CP is also illustrated in Fig. |l|. Ne- 
glecting the possibility of inhomogeneous phases, the black dotted line shows the first order phase 
transition when limiting to homogeneous phases. The shaded domain then depicts where inho- 
mogeneous phases aie energetically preferred. It is enclosed by two second-order phase transition 
lines that meet at the CP, and it covers the first order phase transition line. 

We do not want to address the general question whether an inhomogeneous phase with a higher 
dimensional modulation could become favored in the vicinity of the CP, but it may very well be that 
the one-dimensional modulations are generally preferred close to the CP as numerically confirmed 



in Ref. [23] for the analogous case of inhomogeneous phases in paramagnetic superconductors. 



3. Phase diagrams allowing for inhomogeneous ground states with a 
one-dimensional modulation 

In order to determine the mean-field thermodynamic potential without further approximations, we 



evaluate Eq.(|2.3[) to 

a(r,Ai;M(x)) = _££Trz,,i/Log(^-(/«„+//MF-M)j+:^ Jv 4G^' 

where we introduced the energy-spectrum {En} of the Dirac-type Hamiltonian 

HMF = -if7'di + fM{x). (3.2) 

Limiting to order parameters with a one-dimensional modulation M(x) = M(z), we can make use 
of two connected properties. First, we use the conserved momentum pj^ in the perpendicular x- 
and y-direction combined with Lorentz symmetry and express all eigenvalues through the set {A } 



at Pj^ = (see Ref. |14] for details), giving 



(271)2 \ \ IT j j V Jv 4G, 

+const. . (3.3) 

Second, we observe that the Hamiltonian Hmf for the case pj^ = can be cast into the form 

Hmf-w = (^"'" h,d) ' ^^'"^^ 

where Hmf-w is the Hamiltonian of the ON model for the same order parameter M{z). Since 
the spectral densities for the considered inhomogeneous phases are analytically known, we can 
evaluate the thermodynamic potential and in particular the sum over the eigenvalue spectrum. 
Using the gap-equation it can furthermore be shown that self-consistent solutions, i.e. local minima 
of the thermodynamic potential, in the GN model can generically be used to find self-consistent 
solutions in the NIL model [Il4|l. 
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^^q [MeV] [MeV] 

Figure 2: Left: Structure of the NJL phase diagram in the chiral limit as a function of temperature T and 
quark chemical potential pLq for Mq = 300MeV. The black (short-dashed) line indicates the second order 
phase transition from chirally broken to restored phase, the red (solid) line the first order phase transition 
and the bullet the critical point. The spinodal region is enclosed by the blue (long-dashed) lines. Right: 
Same plot as on the left including the orange (shaded) domain where the energetically preferred ground state 
is inhomogeneous. 




260 280 300 320 340 360 380 400 260 280 300 320 340 360 380 400 

liq [MeV] ^lq [MeV] 



Figure 3: Left: Same plot as on the right of Fig. 2, now including results forM^ — 350MeV (upper branch) 
and Mq — 250MeV (lower branch). Right: Same plot as on the right of Fig. 2, now including the domain of 
inhomogeneous phases for m = 5MeV and m — lOMeV. Branches with critical points at smaller temperature 
T and larger quark chemical potential pLq correspond to larger current quark masses m. 



Consequently, it is possible to use the analytically known solutions in the l-i-l dimensional GN 
model for exploring the NJL phase diagram including inhomogeneous phases with a one-dimensional 



modulation. The form of the order parameter of interest here has the specific form [16] 



M{z) = Afvsn(Z7|v)sn(Az|v)sn(Az + Z;|v) + ''"^^'^,^,'|"f'^^ ^ , (3.5) 

where A is a scale parameter and sn, cn, dn are elliptic Jacobi functions with elliptic modulus \/v. 
Physically, it describes lattices of equidistant domain-wall solitons and it is worth noting that it also 
parametrizes self-consistent solutions at finite current quark masses m. 

As already noted in the context of the GL expansion, we have to introduce a regularization scheme 
in order to make a specific calculation of a phase diagram. Since the usually employed regulariza- 
tion schemes in momentum space rely on a conserved three-momentum for the quasi-particles and 



therefore on homogeneous phases, we apply a proper-time regularization and refer to Ref . [ 14 ] for 
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details. We fix the coupling constant and cutoff scale A by choosing the pion-decay constant 
f:^ = 88MeV for m = and the constituent quark mass in the vacuum Mq = 250, 300, 350MeV. 
Minimizing the thermodynamic potential in the parameters A, v and b for each /i and T , we obtain 
the phase diagrams shown in Fig. ^ and|3[ 

As suggested by the generalized GL expansion we find that the first order phase transition is re- 
placed by two second order phase transitions enclosing a domain of inhomogeneous phases. The 
location of the CP depends on the model parameters and the region of inhomogeneous phases seems 
strongly correlated with the strength of the first order phase transition when limiting to homoge- 
neous phases. To illustrate the latter we also show the extend of the spinodal region, where the 
thermodynamic potential possesses more than one local minimum. The picture stays qualitatively 
the same when including finite cunent quark masses, only the second order phase transition in the 
case of homogeneous phase transition is rendered a cross-over here. 



4. Inhomogeneous phases in the QM model 

The NJL model regularized by a proper-time regularization and adjusted to chiral condensate and 
pion decay constant is known to give constituent quark masses of order 200MeV in the vacuum. 
Hence it gives an undesired phenomenology with regard to the QCD phase diagram, mainly be- 
cause quasi-particles will start forming a Fermi surface at ~ Mq. This is phenomenologically 
unacceptable for /i < (Ma? — B)/3 ~ 308MeV, where is the nucleon mass and B the binding 
energy of nucleons in nuclear matter. In the study of the phase diagram we have therefore chosen 
to fix /;;: to its phenomenological value and varying the value of M^. Since we haven't found a 
regularization that avoids these problems, we also discuss a model that is very similar to the NJL 
model and where the issue of the regularization scheme can be surpassed: The linear sigma model. 



which in this context is usually named QM model [ 24 , 



The Lagrangian of the QM model with Nf = 2 and Nc = 3 is given by 

^QM = v{ifd^,-g{o + iY5r''K"))Y-Uio,K"), 
U{oX) = --{d^od^o + d^K"d^n") + ^{o'^ + K"%" -v^f -co, (4.1) 

1//^ is again the 4A^/A'c -dimensional quark spinor, o the scalar field of the a-meson and n" the 
pseudo-scalar fields of the pion triplet. In mean-field approximation we treat the fields o and as 
classical and replace them by there expectation values [ 24 , 2^] . Furthermore we can use low-energy 



relations to connect the parameters c, g, X and with hadronic observables. We will express those 
by the pion-decay constant the constituent quark mass in the vacuum Mq, the pion mass m;^ 
and a-meson mass via (a) = f^, {n") = 0, c = mlf^, g = Mq/ f^, X = {m%- m\f / {Ifl) and 

For the thermodynamic potential in mean-field approximation we only include the contributions of 
the fermionic fluctuations and approximate 

Q.qm{T,ii;o{x)) = — ^£Trz),/yLog ( - (/w^+^MKgM-At) J +^ J^U{o{x),n"{x)) , 

(4.2) 
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^^q [MeV] [MeV] 

Figure 4: Left: Phase diagram for the QM model in the chiral limit with the red (dashed) line indicating the 
first order phase transition when limiting to homogeneous phases only. Inhomogeneous phases are preferred 
in the orange (shaded) domain enclosed by black (solid) Hnes. Right; S ame plot as on the left, now for 
niji ~ 138MeV and with the bullet showing the critical point. 



where a(x) is taken to be the only non- vanishing expectation value and the Hamiltonian reads 
Hmf,qm = —iy^ydi + y^gG{x). With the identification M(x) = g(j{x) we can therefore evaluate the 
functional trace-logarithm for the same inhomogeneous phases as in the case of the NJL model. In 
general those phases need not to be self-consistent solutions of the QM model, but for the purpose 
of checking whether an inhomogeneous phase is preferred compared to any homogeneous phase, 
this approach is sufficient. 

The QM model is renormalizable which means that the divergences in the functional trace-logarithm 
can in principle be absorbed by the model parameters. Instead of a proper renormalization we will 



however follow Refs. [g4[ g5|], where it has been assumed that the zero temperature contribution 
can well be approximated by y JyU (a(x), 7r"(x)) with the parameters directly adopted to pheno- 
menology. Choosing the model parameters through fji = 93MeV, = 300MeV, = 600MeV 
and considering the chiral limit mj^ = as well as the physical point with pion mass rrijc = 138MeV, 
we obtain the phase diagrams shown in Fig. ^. The model in the applied approximation does not 
have a critical point in the chiral limit and, as maybe expected from the results obtained in the NJL 
model, the first order phase transition from the chirally broken to the restored phase is replaced 
by two second order phase transitions enclosing a domain where inhomogeneous phases are pre- 
ferred. For the physical value of the pion mass, we find a CP and the qualitative structure of the 
phase diagram is the same as in the NJL model. 
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